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Optimal Wavelets for Electrogastrography 
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Abstract 

Matching a wavelet to class of signals can be of interest in feature detection and classification based on wavelet representation. 
The aim of this work is to provide a quantitative approach to the problem of matching a wavelet to electrogastrographic (EGG) 
signals. Visually inspected EGG recordings from sixteen dogs and six volunteers were submitted to wavelet analysis. Approximated 
wavelet-based versions of EGG signals were calculated using Pollen parameterization of 6-tap wavelet filters and wavelet compres¬ 
sion techniques. Wavelet parameterization values that minimize the approximation error of compressed EGG signals were sought 
and considered optimal. The wavelets generated from the optimal parameterization values were remarkably similar to the standard 
Daubechies-3 wavelet. 
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1 Introduction 

Cutaneous recordings of gastric electrical activity (GEA) known 
as electrogastrography (EGG), can play a major role in the di¬ 
agnosis of gastric motility disorders (TJ. Because of its low-cost 
and non-invasiveness, the EGG technique has great appeal as a 
clinical tool and has been related to various gastric motility ab¬ 
normalities 0. Multiple studies have been conducted in other 
to analyze EGG recordings. Although signal processing of EGG 
signals has been considered essential for extracting clinically rel¬ 
evant information, various traditional methods have been utilized 
with limited success 0. 

Recently, advanced signal processing techniques, such as 
wavelets, have been employed to analyze electrogastrograms (4- 
(8). This approach has been used to (i) propose new wavelets that 
can offer a better time-frequency localization of EGG record¬ 
ings HU |5]; (ii) perform noise detection in EGG signals 0; 
(iii) cancel artifacts related to stimulation 0; and (iv) charac¬ 
terize global gastric electrical dysrhythmias 0. An important 
aspect of wavelet analysis is related to designing a wavelet that 
matches a class of signals. Although wavelet matching can be 
of great importance for detection and classification 0. wavelets 
that match EGG signals have not been systematically sought. 
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The present study addresses the problem of finding a wavelet 
that best “matches” the waveshape of EGG signals in basal state. 
Although there are numerous issues concerning the choice of 
wavelet for signal analysis 0 , generally, a wavelet can be re¬ 
garded as best suited to a class of signals if the latter can be 
represented by as few wavelet coefficients as possible EBED. 
Thus, wavelets which resemble the waveshape of the signal un¬ 
der analysis are often selected. 

In the framework of the proposed research methodology, 
an optimal wavelet is sought that can adequately represent a 
wavelet-compressed EGG signal at a given compression ratio. 
The optimality is detected by minimizing an error measure be¬ 
tween the original signal and its compressed version, subject to 
the choice of wavelet. If, for a given wavelet, the error associ¬ 
ated with the compressed signal were minimal, then its wavelet 
coefficients were considered to best represent the original signal. 
Therefore, the selected wavelet would more effectively “match” 
the signal under analysis when compared to other wavelets in 
consideration 02). 

Consequently, the ultimate aim of this study is to quantita¬ 
tively determine a wavelet suitable for the analysis of basal EGG 
recordings in canine and human models. 

2 Methods 

2.1 Experimental Setup 

2.1.1 Canine Experiments 

After a laparotomy and the installment of six pairs of internal 
subserosal stainless steel wire electrodes into the antral gastric 
wall of sixteen acute dogs (seven female and nine male), the 
abdominal wall was closed and five neonatal electrodes (Con- 
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2.1.3 Signal Preprocessing 



Since both canine and human recordings were of significant du¬ 
ration, the raw EGG data were intermittently contaminated with 
a multitude of artifacts, including: (i) motion artifacts; (ii) spon¬ 
taneous variations in electrode potentials; (iii) respiration; (iv) 
signal saturation during recording; (v) electrocardiac activity; 
and fvi) loss of signal during recording. Usually these arti¬ 
facts appeared simultaneously in all recording channels. Some 
of these noisy patterns were visually evident (e.g., iv and vi) and 
could be easily identified and discarded 0. This practice has 
been recommended before in order to obtain a more reliable sig¬ 
nal for subsequent analysis m. 

Therefore, for each subject, a 10-minute time interval of 
channel-synchronized data was manually selected. These data 
were considered to be free from identifiable noise patterns. 

2.2 Signal Analysis 
2.2.1 Wavelet Compression 


Figure 1: Internal (a) and cutaneous (b) electrode positioning in 
canine experiments. Various electrode combinations were used 
for the GEA (c) and the EGG (d) recordings. The electrode com¬ 
bination for the EGG recordings in human experiments was sim¬ 
ilar. 

med, Andover Medical, Haverhill, MA, USA) were placed cuta- 
neously along the projection of the gastric axis. These five elec¬ 
trodes were grouped to form eight EGG channels. In addition, 
the setup provided six internal GEA channels. However, in the 
present study only the EGG channels were processed, while the 
internal GEA recordings were used as a visual reference only to 
verify that normal electrical activity was present. The electrode 
combination set and a diagram of the physical location of the 
electrodes are depicted in Figure [T] 

Thirty-minute EGG recordings were performed in the basal 
state. The captured EGG signals were conditioned by a 0.02- 
0.2 Hz low-pass first order Butterworth active filter. After ampli¬ 
fication, 12-bit analog-to-digital conversion was performed using 
a sampling frequency of 10 Hz and LABMASTER 20009 16- 
channel analog-to-digital converter (Scientific Solutions, Van¬ 
couver, BC, Canada). 

2.1.2 Human Experiments 

Using a similar 8-channel EGG configuration, one-hour record¬ 
ings from six normal volunteers (two female, four male) in post¬ 
prandial state (500 Kcal, 52% carbohydrates, 19% proteins, and 
29% fat) were obtained. The average body mass index for the 
volunteers was 22.2 kg-m -2 (SD 3.0 kg-m~ 2 ). Signal condition¬ 
ing, amplification and digitization process similar to the ones uti¬ 
lized in the canine experiments were implemented. 

All experiments were approved by the Animal Welfare Com¬ 
mittee and the Ethics Committee at the Faculty of Medicine, Uni¬ 
versity of Alberta. 


In a discrete-time formalism wavelet transforms are performed 
via the Fast Wavelet Transform using Mallat’s pyramid algo¬ 
rithm for decomposition (forward transform) and reconstruction 
(inverse transform) m. 

Let x be a discrete signal with N = 2 J points (a sampled ver¬ 
sion of the analog signal x(t)). The discrete wavelet transform 
(DWT) of x is computed in a recursive cascade structure consist¬ 
ing of decimators [2 and complementing filters h (low-pass) and 
g (high-pass), which are uniquely associated with a wavelet 1H. 
Fig. [2] depicts a diagram of the filter bank structure. 

At the end of the algorithm computation, a set of vectors is 
obtained {d i, cU,..., dj,... , d ,/ 0 , a / 0 } , where Jo is the number of 
decomposition scales of the DWT. This set of approximation and 
detail vectors represents the DWT of the original signal. Vec¬ 
tors d, contain the DWT detail coefficients of the signal in each 
scale j. As j varies from 1 to Jo, a finer or coarser detail co¬ 
efficient vector is obtained. On the other hand, the vector a ./ 0 
contains the approximation coefficients of the signal at scale Jo¬ 
lt should be noted that this recursive procedure can be iterated J 
times at most. Usually, the procedure is iterated Jo < J times. 
Depending on the choice of Jo, a different set of coefficients can 
be obtained. Observe that the discrete signal x and its DWT have 
the same length N. The inverse transform can be performed us¬ 
ing a similar recursive approach flOll . Generally, a signal can be 
subject to various wavelet decompositions. The analysis depends 
on (i) the choice of wavelet (filters h and g); and (ii) the number 
of decomposition levels (scales) Jo. 

A wavelet-based compression scheme aims to satisfactorily 
represent an original discrete signal x with as few DWT coeffi¬ 
cients as possible mDEna. One simple and effective way of 
doing that is to discard the coefficients that, under certain crite¬ 
ria, are considered insignificant. Consequently, the signal recon¬ 
struction is based on a reduced set of coefficients muni. 

In the present work the classic scheme for non-linear com¬ 
pression was used QU. This procedure considers an a posteriori 
adaptive set, which keeps M wavelet transform coefficients that 
have the largest absolute values. A hard thresholding was used 
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Figure 2: Wavelet analysis filter bank. The signal is iteratively 
decomposed through a filter bank to obtain its discrete wavelet 
transform. 


to set the remaining coefficients to zero. The number of coeffi¬ 
cients M to be retained was determined according to the desired 
compression ratio CR, which was defined by 


CR = —, 
M 


(1) 


where N and M are the number of wavelet transform coefficients 
of the original and the compressed signals, respectively. 


2.2.2 Measurement of Distortion 


To further the analysis, it is necessary to introduce an error mea¬ 
sure to compare the original discrete signal x with its reconstruc¬ 
tion x. Several measures that allow the evaluation of the effect 
of compression schemes have been suggested lfl8ll . However, 
one of the most commonly used is the Percent Root-mean-square 
Difference (PRD) IT5llT6l[T8l. which was utilized in the present 
study as a measure of distortion in the compression scheme. The 
PRD of two signals, x and x, both of length N, is defined by: 


PRD(x,x) 



x 100%. 


( 2 ) 


2.3 Choice of Parameters 
2.3.1 Number of Scales 

In order to select the number of scales Jq £ {1,...,J} of the 
wavelet transform decomposition, the following criterion was in¬ 
troduced: Jo was chosen so that the coarsest approximation scale 
had a pseudo-frequency close to the EGG dominant frequency f c 
of 4-6 cycles per minute for the canine subjects fl9l and 3 cycles 
per minute for the humans. 

The pseudo-frequency / pse udo of a given scale j is 

/pseudo = -Ajr, 7 = 1,2,...,/, (3) 

J - T s 

where T s is the sampling period (0.1 s) and fy is the center fre¬ 
quency of a wavelet (the frequency that maximizes the mag¬ 
nitude of the Fourier transform of the wavelet) (331 ■ Conse¬ 
quently, a scale Jq was selected which minimized the difference 
(./pseudo — fc )- TableQ]shows the number of decomposition levels 
for some common wavelets. 


Table 1: Number of decomposition scales Jq for some wavelets 


Wavelet 


Jo 

Canine 

Human 

Haar 

7 

8 

Daubechies-2 

6 

7 

Daubechies-3 

7 

7 

Coiflet-1 

7 

7 


2.3.2 Compression Ratio 

A compression ratio set CR £ {3,5,7,10} was selected and a 
matching procedure was carried out aiming at optimizing the 
choice of wavelet. 

2.4 Optimization of the Wavelet Choice 

In the context of the present study, a wavelet was sought that 
minimized the PRD between the original EGG signal and its re¬ 
construction for a given compression ratio. 

However, the abundance of wavelets m makes such ap¬ 
proach prohibitive. As a result, some constraints on the choice 
of wavelet were introduced. 

It is well known that wavelets can be generated from discrete 
finite impulse response (FIR) filters m. In the present work, 
the analysis was limited to wavelets generated by FIR filters with 
length no greater than six coefficients. In this subset of wavelets 
one may find Haar, Daubechies-2, Daubechies-3, and Coiflet-1 
wavelets, to name the most popular ones m. 

This restriction is quite convenient, since all FIR filters of 
length up to six that can be utilized to generate wavelets have 
simple parameterizations of their coefficients |[2()1 . For exam¬ 
ple, Pollen parameterization of 6-tap wavelet filters ED has two 
independent variables (a,b) £ [—n,n] x [—n,n\. Varying these 
two parameters, a filter that generates a new wavelet can be de¬ 
fined. Consequently, the Pollen parameterization defines a plane 
on which every point is connected to a wavelet |2TI . 

Using the discussed compression scheme, one can compute 
a PRD value for each wavelet generated from a point with coor¬ 
dinates (a,b) on the parameterization plane. Doing so, a surface 
can be defined by the points (a,b, PRD). Thus, the minima of 
this surface correspond to the point coordinates ( a,b) that gen¬ 
erate a wavelet with good “matching” properties, since the PRD 
values at these minima are small. 

As a result, a set of point coordinates ( a, } bi ) could be deter¬ 
mined on the parameterization plane, which minimizes the PRD 
for each EGG recording i. Fig. [3 shows typical surfaces gener¬ 
ated for basal canine and human EGG signals. 

Taking the mean value of the minima, the best wavelet pa¬ 
rameterization (a*,b*) could be defined. Thus, (a*,b*) gener¬ 
ates a wavelet that on average “matches” best the normal EGG 
recordings. 

3 Results 

The optimal values of the wavelet parameterization were deter¬ 
mined for the selected compression ratios (Table [2j. The val¬ 
ues for compression ratio of 3 were associated with the wavelets 
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a a 

(a) CR = 3 (b) CR = 10 



a a 

(c) CR = 3 (d) CR = 10 


Figure 3: Plots generated after computing the PRD surface for all possible wavelets on the parameterization plane, using a canine 
EGG signal with CR = 3 (a), and CR = 10 (b). A representative human EGG signal was used to build PRD surfaces with CR = 3 (c), 
and CR = 10 fd). The minimum value is depicted by a circle (o). The coordinate points that correspond to Haar (*), Daubechies-2 
(x), Daubechies-3 (•), and Coiflet-1 (A) wavelets are shown. The axes are normalized by n. 
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Table 2: Optimal values of the parameterization 


Canine^ Human^ 



a* 

b* 

a* 

b* 

3 

0.4329 

—0.2608 

0.3976 

-0.2335 

5 

0.4323 

-0.2638 

0.4293 

-0.2550 

7 

0.4323 

-0.2700 

0.4276 

-0.2551 

10 

0.4293 

-0.2736 

0.4279 

-0.2600 


“Based on 16 canine subjects. 
Glased on 6 volunteers. 



Figure 4: Comparison between the obtained optimal wavelets for 
canine (dotted curve) and human (solid bold curve) EGG signals 
and the standard Daubechies-3 wavelet (solid thin curve) using 
the optimal parameterization for compression ratio of 3. The 
similarity between Daubechies-3 and the optimal EGG wavelets 
is clearly evident. 

depicted in Fig. 0] It is worth mentioning that the proposed 
wavelets were very similar to the classic Daubechies-3 wavelet. 
The correlation coefficient between Daubechies-3 and the opti¬ 
mal EGG wavelet was 0.996 and 0.965, for canine and human 
signals respectively. Since Daubechies-3 is (i) very similar to 
the proposed wavelet; and (ii) easily available in many software 
packages (e.g., Matlab (The Mathworks, Inc., Natick, MA, 
USA)), it can be chosen instead of the proposed wavelets. Previ¬ 
ous empirical findings Q confirm this observation. 

4 Conclusion 

The problem of finding optimal wavelets to “match” EGG sig¬ 
nals was quantitatively addressed. The proposed wavelets can be 
considered as tools to further EGG signal analysis. Moreover, 
the suggested methodology opens an avenue towards the classi¬ 
fication of electrogastrograms based on the PRD value of their 
wavelet compressed version, either applying the obtained opti¬ 
mal wavelets or the standard Daubechies-3. 
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